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I. INTRODUCTION 



We formulate the nonlocal NJL model with a momentum dependent 
constituent quark mass and calculate pion light cone wave functions of 
twist 2 and 3. The leading twist wave function is not asymptotic and 
_Cj ' agrees well with the new CLEO data. Normalization conditions for 

the twist 3 wave functions are used to calculate the quark condensate. 
A prescription to calculate the gluon condensate is proposed. The 
If} ' numerical value of the gluon condensate nicely agrees with the phe- 

£N) , nomenological value, whereas the quark condensate is larger than the 

phenomenological value of —(250 MeV) 3 . The relation between the 
fcy moments and mixed condensates are used to estimate the mixed 
. quark-gluon condensate of dimension 5. 
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In this short note we shall describe a simple and tractable model for the pion light 
cone wave functions which is based on the instanton model of the QCD vacuum. 
Hadron light cone wave functions were theoretically introduced more than 20 years 
Qh! &E° 0- Pi- Recently the analysis of Ref. || based on the latest CLEO measurements 

0| put some limits on the expansion coefficients of the axial- vector (AV) pion wave 
" function in terms of the Gegenbauer polynomials. This analysis indicates that 

the pion wave function measured at Q 2 = 1.5 — 9.2 GeV 2 is neither asymptotic 
(j>^y (u) = 6 u(l — u) (with u being the fraction of the pion momentum carried 
by the quark) nor of the form proposed by Chernyak and Zhitnitsky in 1977 
■ ficzi 11 ) = 30 u(l — u)(l — 2m) 2 . These two wave functions together with a typical 

prediction of the present model are shown in Fig. la. In Fig. lb we show the 95% 
and 68% confidence level contour plots in the — &4 parameter space from the 
analysis of Schmedding and Yakovlev (Fig. 6 in Ref. ||) together with the values of 
d2 and ci4 for , (pQ^ and various parameters of the present model. 

The instantanton model, after integrating out gluons and performing the 
bosonization, reduces to a simple Nambu-Jona Lasinio type model where the quarks 
interact nonlocally with an external meson field U P,JL0| : 



-^^k)^Wk)u^{k i)Vmjm (i) 

and U 15 can be expanded in terms of the pion fields: 

U^=l + ^--y 5 r\ A -^\ A + ... (2) 
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Here = 93 MeV and M(k) = MF 2 (k) is a momentum dependent constituent 
quark mass which also plays a role of the pion-quark coupling. Let us note that in 
the instanton model both quark and gluon condensation occur at the same scale (1q 
which is associated with the average instanton size 1/p = 600 MeV. 

In principle F(k) has been calculated in the instanton model in the Euclidean 
space time. Here, following Refs. |y],[^), we will perform calculations directly in 
the Minkowski space. To this end we shall choose a simple pole formula 

/ A 2 \ n 

m= - to A2| ... (3) 



fc 2 -A 2 

which for n ~ 2 — 3 and for k 2 < reproduces the k dependence obtained from the 
instantons reasonably well |L2| (see Fig. 2a). Here M = M (0) is a model parameter 
which we choose to be of the order of 350 MeV. 
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FIG. 1. Left: Asymptotic and Chernyak-Zhytnitsky leading twist pion wave functions 
together with a typical wave function from the present model. Right: The parameter 
space (02,04) of Ref. [6]. Black dots represent different model predictions, solid contour 
corresponds to 68% confidence level, whereas the dashed one to 95% 



As we shall see, the model is technically very simple and allows to calculate 
pion wave functions (not only the axial-vector, but also the pseudo-scalar (PS) and 
the pseudo-tensor (PT) ones) analytically up to a numerical solution of a certain 
algebraical equation of the order An + 1. Given this simplicity it is of importance 
to perform various tests in order to gain confidence in the model as well as to find 
its limitations. In this paper we provide 4 kinds of tests. 

First we calculate the leading twist pion wave function and compare with the 
existing data. Next we calculate the non- leading twist wave functions, which are 
normalized to the quark condensate. This allows us to calculate (qq). 

It is important to note that in our approach we calculate not only the u depen- 
dence but also the dependence on the transverse momentum fc^: 
00 1 
4>tt(u) = J dfcy Vtt(w, 0tt(A;t) = J duipniu,^). (4) 



By calculating k\ moments we get the mixed condensate of dimension 5. 

Another advantage of our method is that the analytical expression for the quark 
condensate is given in terms of a Minkowskian integral which in a limit of a constant 
M{k) and k 2 — > —k 2 E reduces to the well known Euclidean form. By comparing 
the two expressions one can by inspection guess a continuation prescription which 
allows to rewrite certain Euclidean integrals as the Minkowskian ones. We use this 
in some respect ad hoc prescription to calculate the gluon condensate (a/ir GG), 
which provides another test of our approach. 
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II. PION WAVE FUNCTIONS IN THE NONLOCAL QUARK MODEL 



We shall be dealing with the leading twist axial- vector (AV), twist 3 pseudo-scalar 
(PS) and pseudo-tensor (PT) wave functions defined as follows Jl3|,p|: 



# V («) 



1 



dr 

iVZFir J n 

— oo 



^_ e -ir(2u-l)(nP) ( Q | ^{n T )^ l5 ^ (_ nr ) |7r+(p)) , 



(5) 



(h 



- t , -e— ( 2 "- 1 )(" p )(0|^(nr) l75 ^(-nr)|^+(P)), 
\/2{qq) .1 ir 



u OO 

_^[ dw f ^e--( 2 - 1 )^n«P^0|Vi(nr)^75V'(-nr)|7r+(P)>. 
V2 (to) J J it 

— oo 



where we have chosen n = (1,0,0,-1) as a light-cone vector parallel to z M = rn^ 
and n = (1,0,0,1) parallel to P M . All three wave functions are normalized to 1. 
The normalization condition for <p^ s and yield (j>^ s therefore the expression for 
(qq), whereas normalization of $f is used to fix the model parameter A for given 
M and n. 

Technically speaking all three wave functions (||) arc given in terms of a loop 
integral with a momentum dependent quark mass M(k), which also acts as a quark- 
pion coupling. In order to calculate the loop integral we have to find zeros of the 
propagators which are generically of the form 



fc 2 - M 



A 2 



fc 2 - A 2 



4)1 



ie = 0. 



Equation (^|) can be conveniently rewritten as: 



(6) 



(7) 



where z — k 2 /A 2 — 1 + ie and fx 2 
d 4 k = l/2dk+dk_-d 2 k T where fc' 1 



- M 2 /A 2 . In the light cone parametrization 
k+h^/2 + k-n^/2 + fc£. Since k+ = uP+ is 
fixed, equation (^|) should be understood as an equation for k~ . Generally, equation 
(0) has 4n + l complex solutions which in the following will be denoted as Zi. These 
solutions depend on the specific value of /z 2 and have to be calculated numerically. 

Here one faces immediately the problem how to choose the integration contour 
in the complex fc~ plane. The prescription is very simple and has been at length 
discussed in Ref. |l2[ . As a result the dk~ integrals yield real wave functions which 
vanish for u outside the region < u < 1. Moreover for A — > oo, i.e. for a 
constant M (fc), this prescription reduces in a continuous way to the standard one 
of Feynman. 

With this prescription the calculations are rather straightforward and we obtain: 



N c M 



in+l 

(2tt) 2 A 2 F 2 ^ 



fifk 



+ (1 - u)z? 



t + 1 + UZi + (1 — u)Zk ' 



(8) 



^ S (u,k 2 T ) 



N c M 
(2tt) 2 (qq) 



4n+l 



J2 M*'- 



v 3n „3n 



(i 



•- + Zfc ' 



m 2 4 



i,k=l 



t + 1 + UZi + (1 - u)z k 



(9) 
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oat ]ur\2 4n+l 

^r(«,^) = / : ° a ^\ E /iA^ B ^Mi+t+u* + (i-u)* fc ). (io) 
( 2?r ) w; i,fc=i 

(where t = fc^/A 2 ). Factors fi obey the following properties: 

4n+l 4n+l ( for to < 4n 

/< = n — > E = < (ii) 

k=i z% Zk i=i 1 for to = 4/i 

fc^j v 

which are crucial for the convergence of the dt integrals. 

It is now straightforward to perform either the dt integration in order to get cf> w , 
or the du integration to get the kx -dependent functions <p w . 



III. PROPERTIES OF THE PION WAVE FUNCTIONS 
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FIG. 2. Left: F(k) for Euclidean momentum k 2 < 0, for n — 1 (dashed), 5 
(dashed-dotted) and for the instanton model (solid). Right: Axial- vector pion wave func- 
tion for M = 350 MeV and for n = 1 (dashed) and 5 (dashed-dotted) together with the 
asymptotic one (solid) 

In order to study the model dependence on the choice of M and n we have 
calculated pion wave functions for M = 325 — 400 MeV and n = 1 — 5. The cutoff 
parameter A was adjusted by imposing the normalization condition on 4^ V • I n 
fact, as discussed in Ref. |Q, the leading twist pion wave function 4>£ v (u) does not 
change any more if we increase n above 5. On the other hand for n > 5 the cutoff 
function (||), if continued to the Euclidean metric, starts to deviate significantly 
from the one obtained in the instanton model. Therefore we have chosen to work 
with n max = 5. In Figs. 2b and 3 we have plotted (f)^ v , 4>^ s and <j)^ T for M = 350 
MeV and n = 1,5. 
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FIG. 3. Pseudo scalar (left panel) and pseudo-tensor (right panel) pion wave function for 
M = 350 MeV and for n — 1 (dashed) and 5 (dashed-dotted) together with the asymptotic 
one (solid) 



Let us shortly summarize our findings. The axial- vector wave function, 4>^ v , 
vanishes at the end points as u n (or (1 — u) n ) and shows a plateau around u = 0.5 
with a small dip for n = 5. It differs from the asymptotic wave function and, 
as seen from Fig. lb, the best agreement with the recent analysis of the CLEO data 
is obtained for M = 325 MeV and n = 2 - 5 or M = 350 MeV and n = 2. The 
fact that the true pion distribution amplitude may be broader than the asymptotic 
one has been already pointed out in Ref. p5| . Such a behavior was found then in 
Ref. [|l6| where not only the nonlocality (within the sum rules approach) but also 
the radiative corrections have been taken into account. 

The pseudo-scalar pion wave function, <fi^ s , was calculated within the QCD sum 
rules in Refs. 13 14|]. It had a u— shape and did not vanish at the end points. In 



our case (j)^ s vanishes at the end points for n > 1. Indeed 



y in in , 1 3n+l in , 1 in in+l 2 n n 

-i z k t 2 Z i Z k t 2 i k " * k 

(12) 



(and similarly for u = 0) which is equal to due to the property (|ll| ) except for 
n = 1 where in + 1 = An. Interestingly, the vanishing of (f>^ s for u = 0, 1 is 
correlated with the nonconvexity of 4>^ v at the end points, which, as stated above, 
behaves like u n (or (1 — it)") for it — * (or 1). In any case 0^ >s differs from its 
asymptotic form (f>^ s s = 1. 

Both pseudo-scalar and pseudo-tensor wave functions show stronger rt dependence 
than 4>n V ■ F° r ^ = 1 <^ T coincides with the asymptotic expression = (j)^ , 
while for n = 5 its is depleted at the end points and peaked in the center. 



IV. CONDENSATES 



Since the model parameters are fixed by the normalization of the axial-vector 
wave function we could use the normalization condition for (p^ s ov (j)^ T to calculate 
the quark condensate. The results are presented in Table 1. We see that the quark 
condensates obtained from the two normalization conditions do not coincide. In 
fact, for almost all model parameters considered, we find 

ff(qq)ps/(M)pT^^- (13) 



5 



In absolute values the quark condensate calculated within our model overshoots the 
phenomenological value of —(250 MeV) 3 . This is mostly due to the rather poor 
convergence of the dk\ integrals (^) . Indeed for large k\ : 

(-. \ 4n+l / , \ In / . \ 2tl 

This is also the reason of rather strong n dependence of ^^and <p^ T . 

The Euclidean formula for the gluon condensate in the instanton model of the 
QCD vacuum reads ||: 

a \ f d 4 k E M 2 (k E ) , s 

-GG) = 32V C / —^—^ ; r9 , . . 15 

7T / C J 2tt 4 k% + M 2 (k E ) 



Apart from the numerical factor in front it differs from (qq) by an additional power 
of M(k E ) in the numerator. In Ref. |T^] we have suggested the continuation pre- 
scription of (|l5|) to the Minkowski metric with the result 



7T / (2tt) 2 7 fr' i + 1 + + (1 - u)z k v ; 

Numerical result (Table 1) depends very weakly on n and is compatible with the 
phenomenological value |j: (a/ it GG) = (393±§| MeV) 4 . 





TABLE I. 


Condensates for M 


= 350 MeV 




n A 


(fGG) 


(QQ) ps 


(q<i)pt 


(igqcr- Gq) AV 


1 1156 MeV 


(399 MeV) 4 


-(318 MeV) 3 


-(357 MeV) 3 


-(553 MeV) 5 


5 2819 MeV 


(389 MeV) 4 


-(271 MeV) 3 


-(301 MeV) 3 


-(475 MeV) 5 



Soft pion theorems provide link between dynamical objects like the light cone 
wave functions [p]-pj and static properties of the physical vacuum |@. It has 
been shown in Refs. ^J,^) that moments of ^(k^) are given in terms of the mixed 
quark-gluon condensates 



/,2\ _ 5 (igqcr ■ G q) , 2 \ _ 1 (igqcr ■ G q) 
\ t) av ~ 36 (^) ' Wps - 4 ^ ■ ^7) 



Here G£„ is a gluon field strength and a ■ G = a^G^, G MV = A Q /2 G£„. 

Unfortunately the ratio (jt^). v /(k? r ) ps ~ 5/9 which follow from ( p7| ) is not 
reproduced within our approachn due to the slow convergence of the dk\, integration 
in the case of (Jc^) ps . In order to estimate the value of the mixed condensate 
of dimension 5, (igqcr -Gq), we choose therefore the first equation of Eqs.(|l7|). 
Interestingly, for the parameters which are closest to the original instanton model, 
M = 350 MeV and n = 2, we get —(493 MeV) 5 in perfect agreement with the direct 
calculation of (igqcr ■ Gq) in the instanton model jl9| which gives —(490 MeV) 5 . 



x It is of the order of 0.2 instead of 0.54 



6 



V. SUMMARY AND OUTLOOK 



The nonlocal NJL model with the momentum dependent constituent quark mass 
has been applied to calculate pion light cone wave functions fll| . It gives a satis- 
factory description of the leading twist AV wave function, whereas for the twist 3 
wave functions we find a somewhat larger sensitivity to the model parameters. 

Present prescription can be easily extended to describe kaon wave functions with 
an explicit symmetry breaking due to the non zero current strange quark mass. 
Also two meson generalized parton distributions both for pions and kaons can be 
easily calculated. By crossing symmetry one can also apply our method to calculate 
the skewed distributions and structure functions |2^ , |23[| . 

On the theoretical side one has to investigate more closely the PCAC relation 
within the present approach. It is known that the properly defined currents should 
include additional terms with respect to those considered here po|j2l| . Although 
these new terms are not unique and suppressed by the instanton packing fraction, 
their influence on our results should be investigated. 
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ing. This work was partially supported by the Polish KBN Grant PB 2 P03B 019 17. 
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